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Abstract 

Einstein- Yang-Mills-dilaton theory possesses sequences of neutral static spher- 
ically symmetric black hole solutions. The solutions depend on the dilaton cou- 
pling constant 7 and on the horizon. The SU(2) solutions are labelled by the 
number of nodes n of the single gauge field function, whereas the SO (3) solu- 
tions are labelled by the nodes (ni,n2) of both gauge field functions. The S0(3) 
solutions form sequences characterized by the node structure (j, j -|- n), where j 
is fixed. The sequences of magnetically neutral solutions tend to magnetically 
charged limiting solutions. For finite j the S0(3) sequences tend to magnetically 
charged Einstein- Yang-Mills-dilaton solutions with j nodes and charge P = \/3- 
For j = and j ^ 00 the S0(3) sequences tend to Einstein-Maxwell-dilaton 
solutions with magnetic charges P = -v/3 and P = 2, respectively. The latter also 
represent the scaled limiting solutions of the SU(2) sequence. The convergence of 
the global properties of the black hole solutions, such as mass, dilaton charge and 
Hawking temperature, is exponential. The degree of convergence of the matter 
and metric functions of the black hole solutions is related to the relative location 
of the horizon to the nodes of the corresponding regular solutions. 
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1 Introduction 



Motivated by higher- dimensional unified theories, such as Kaluza-Klein theory or su- 
perstring theory, where a scalar dilaton field arises naturally, static black hole solutions 
have been found in Einstein-Maxwell-dilaton (EMD) theory [|l], 0. In contrast to the 
static charged black holes of Einstein-Maxwell theory, referred to as Reissner-Nordstr0m 
(RN) black holes, the charged EMD black holes exist for arbitrarily small event horizon, 

> 0. The "extremal" solutions of EMD theory, obtained in the limit xh = 0, possess 
a naked singularity at the origin. 

Recently static black hole solutions have also been studied in Einstein- Yang-Mills- 
dilaton (EYMD) theory 0, |, §, §, 0. SU(2) EYMD theory possesses a sequence of 
magnetically neutral static spherically symmetric black hole solutions, labelled by the 
number of nodes n of the single gauge field function. The solutions exist for arbitrary 
event horizon xh > p, H, ||, H, 0. Here, in the limit xh — regular particle- like 
solutions are obtained P, |, |], |], H, 0. 

For finite dilaton coupling constant 7, the sequence of regular neutral SU(2) EYMD 
solutions tends to the "extremal" EMD solution with the same dilaton coupling constant 
7 and with charge P = 1 for large n H]. For finite event horizon xh and dilaton 
coupling constant 7, the corresponding sequence of neutral black hole solutions of SU(2) 
EYMD theory also tends to a limiting solution for large n. This limiting solution is the 
EMD black hole solution with the same event horizon xh, the same dilaton coupling 
constant 7, and with magnetic charge -P = 1 [0- 

The magnetically neutral black hole solutions of EYMD theory have many features 
m common with those of Einstein- Yang-Mills (EYM) theory |, |10|, |ll], |12|. In fact, 
in the limit of vanishing dilaton coupling constant 7, the EYMD black hole solutions 
approach those of EYM theory. In SU(2) EYM theory, the limiting solutions are the 
RN black holes with charge P = 1 and event horizon xh, for xh > 1 [0, II^ S 



For sequences of black holes with xh < 1 and for the regular sequence IJI^, the limiting 
solutions are different 0, [T|, |r|. 

The non-abelian SU(2) black hole solutions do not possess a global YM charge. Since 
they are characterized not only by their mass but in addition by an interger n, they 
represent (unstable |T8|, |T9|, |20|, |21|, |2^) counterexamples to the "no-hair conjecture" . In 
contrast, SU(2) black holes with non-zero YM charge are uniquely characterized by their 
mass and charge. Indeed, for SU(2) all static spherically symmetric EYM black hole 



solutions with non-zero YM charge are embedded Reissner-Nordstr0m solutions p3 



2^ . This "non-abelian baldness theorem" no longer holds for SU(3), where black hole 



solutions with both Coulomb-like and essentially non-abelian gauge field configurations 
exist p5| , P^ ]. While these black hole solutions correspond to SU(2)xU(l) 



solutions, the properties of genuine SU(3) black hole solutions have remained largely 
unknown. 
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The genuine magnetically neutral static spherically symmetric black hole solutions 
of SU(3) EYM and EYMD theory are based on the SO (3) embedding. These solutions 
can be labelled by the nodes (^1,77,2) of both gauge field functions ^7\, [1^, 0. We 
have constructed SU(3) EYM and EYMD black hole solutions with a small number of 
nodes previously fl^, 0] , and we have identified two sequences of solutions together with 
their limiting solutions 0. These sequences are the sequence (0, n), containing the two 
lowest genuine SU(3) solutions, and the sequence (n, n) of scaled SU(2) solutions. The 
limiting solutions of these two sequences are again EMD solutions with the same event 
horizon and dilaton coupling constant, but with magnetic charges P = a/3 and P = 2, 
respectively. The limiting solutions for 7 = and > -P are the corresponding RN 
black holes. 

For SU(3) EYMD theory it is not a priori clear, in which way the black hole solutions 
with general node structure (ni, 712) assemble to form sequences, converging to limiting 
solutions. One may expect though, that the limiting solutions of the general sequences 
of genuine SU(3) solutions are again magnetically charged solutions. 

Here we construct genuine SU(3) black hole solutions with general node structure 
(ni, n2). We show, that these solutions fall into sequences with node structure (j, j+n), 
and we identify the limiting solutions of the new sequences with fixed finite j with 
the known magnetically charged SU(3) black hole solutions [^, We study the 
convergence of the matter and metric functions and of the global properties for these 
sequences of solutions. 

In section 2 we derive the SU(3) EYMD equations of motion based on spherically 
symmetric ansatze for the metric and the matter fields. We discuss the boundary 
conditions and obtain relations between the metric and the dilaton field. In section 3 
we consider SU(2) EYMD black hole solutions. We investigate their dependence on the 
parameters xh and 7. We relate the degree of convergence of the matter and metric 
functions of the black hole solutions to the relative location of the horizon to the nodes 
of the corresponding regular solutions. We show, that the global properties of the 
solutions, such as mass, dilaton charge and Hawking temperature 0, l^l) converge 



exponentially to those of the limiting solutions. In section 4 we discuss SU(3) EYMD 
black hole solutions. We identify the general set of sequences of genuine SU(3) solutions. 
We demonstrate that the convergence properties of these sequences are analogous to 
those of the SU(2) sequence. In section 5 we present our conclusions. In appendices A 
and B we review the limiting magnetically charged black hole solutions of EMD and 
EYMD theory. We demonstrate that the "extremal" EMD and EYMD solutions agree 
in many respects. 



3 



2 SU(3) Einstein- Yang-Mills-Dilaton Equations of 
Motion 



2.1 SU(3) Einstein- Yang-Mills-Dilaton action 

We consider the SU(3) Einstein- Yang-Mills-dilaton action 

S = Sg + Sm = J Lcv^d'^x + J LuV^d^^x (1) 

with 



and matter Lagrangian 

Lm = -\d,<^d^<^ - e^^^h:i{F,,Fn , (3) 

where 

F^,u = d^A^ - d^A^ - ig[A^, A^] , (4) 
A, = Ix'^A; , (5) 

and g and k are the gauge and dilaton couphng constants, respectively. 

Variation of the action eq. ([^) with respect to the metric g^'^ leads to the Einstein 
equations 

G/iiy = Rfiu — -^QfiuR = SnGTfj,^ (6) 

with stress-energy tensor 



T^v — g^vLM ~ 2 



dg 



= d,^d,<^ - + 2e'^''TT{F,^F,pg"^ - -^g^.F^^F^") , (7) 

and variation with respect to the gauge field and the dilaton field $ leads to the 
matter field equations. 

2.2 Static spherically symmetric ansatze 

To construct static spherically symmetric black hole solutions we employ Schwarzschild- 
like coordinates and adopt the spherically symmetric metric 

ds"^ = g^^dx^'dx'' = -A^Nde + A^-Mr^ + r'^^dO'^ + sin^ Odcj)'^) , (8) 
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with 



N = 1 



2m 
r 



(9) 



For SU(2) EYMD theory, the static spherically symmetric ansatz for the gauge field 
with vanishing time component is, 



^0 

A, 



, 

1 — w{r) 

2rg 



[Cr X r 



(10) 



with the SU(2) Pauli matrices r = {ti,T2,T3). For SU(3) EYMD theory, generalized 
spherical symmetry for the gauge field is realized by embedding the SU(2) or the SO (3) 
generators in SU(3). In the SU(2)-embedding, the ansatz is given by (|10D with r 
replaced by A = (Ai, A2, A3). In the S0(3)-embedding, the corresponding ansatz for the 
gauge field with vanishing time component is 



Ai 



0, 

2 - K{r) 
2rg 



E{r) 
2rg 



{cr X A)i,e; ■ A) 



(11) 



where [ , ]+ denotes the anticommutator, and A = (A7, — A5, A2). For the dilaton field 
we take $ = <^{r). 



2.3 Field equations 

For the above ansatze we now evaluate the tt and rr components of the Einstein 
equations. The metric (Bf) yields the components of the Einstein tensor 



Gtt 



2m'A^N 



SttGTu = -SttGA^NLm , 



(12) 



and 



Gtt 2 A' 
■^2^2 ^ r A 



SnGTr., 



(13) 



where the prime indicates the derivative with respect to r. 

The spherically symmetric ansatze for the fields, eqs. (|10D and ([TT|), yield the tt 
and rr components of the stress-energy tensor, Ttt = —A'^NLm, with 



2 



2k<I> 



Nw'' + — (w'-l 
2r 



(14) 



and T^. 



T 



N' 
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2e 



2k$ 



■W' 
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(15) 
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for the SU(2) embedding, and with 



2k<!> r 



2 g'^r'^ 



and 



Trr = i-LM + + ^{K" + H'') 

for the SO (3) embedding. 

Let us now introduce the dimensionless coordinate x, 



(16) 



(17) 



X 
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(the prime indicating now the derivative with respect to x), the dimensionless mass 
function fi{x) 

(19) 
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m = — ^=m 



the dimensionless dilaton field 0, 

and the dimensionless coupling constant 7, 



7 = k/ vAttG . 



(20) 



(21) 



The choice 7 = 1 corresponds to string theory, whereas 4 + n dimensional Kaluza-Klein 
theory has 7^ = (2 + n) /n |]l| . 

With these definitions we then obtain the dimensionless EYMD equations. For 
SU(2) IP, 1^, H, P] the equations for the metric functions are 



X 



A 



and the equations for the matter field functions are 



{e^^'t'ANw'y = {w^ - l) , 



(22) 
(23) 

(24) 
(25) 



For SO (3) 1^, the equations for the metric functions are 

1 

X 

2 

X 



1 



A' 



1 



T'2\ 



A 



and the equations for the matter field functions are 



e^'"f>ANK'y = ^AK (k^ + - 4 



ie^^^ANH' 



4x2 



AH (h"^ + 7/s:2 - 4 



(AATx^')' = 27^6 



/2 



/2\ 



5X^ 



(26) 
(27) 

(28) 

(29) 
(30) 



With help of eq. (^) for SU(2) or eq. (|27|) for SO (3) the metric function A can 
be eliminated from the matter field equations. Note that the SU(2) equations are 
symmetric with respect to the transformation w{x) — > —w{x), and the S0(3) equations 
are symmetric with respect to both an interchange of the functions K{x) and H{x), 
and to the transformations K{x) —>■ —K{x), and H{x) —>■ —H{x), thereby yielding 
degenerate solutions. 



2.4 Boundary conditions 

Let us now consider the boundary conditions for the black hole solutions of EYMD 
theory. Requiring asymptotically flat solutions implies that the metric functions A and 
fi both must approach a constant at infinity. We here adopt 



A{oo) 



(31) 



thus fixing the time coordinate. For magnetically neutral solutions the gauge field 
functions approach a vacuum configuration, i. e. for SU(2) 



wioo) 



±1 



and for SO (3) 



K{oo) = ±2 , H{oo) = 
K{oo) = , H{oo) = ±2 
For the dilaton field we choose ||, |, H, ^ 

(h(oo) = . 



(32) 

(33) 
(34) 

(35) 
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Magnetically charged solutions require different boundary conditions for the gauge field 
functions. 

The existence of a regular event horizon at xh requires for the metric functions 



2 ' 



and A(xh) < oo. The matter functions must satisfy at the horizon xh 



76 



for SU(2) and 



N'K' 
N'H' 
N'6' 
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Ax^ 



(36) 

(37) 
(38) 

(39) 
(40) 
(41) 



for SO (3). 

On considering the regular solutions, we observe that they satisfy the same boundary 
conditions at infinity as the black hole solutions. However, at the origin, regularity of 
the solutions requires 

/i(0) = , (42) 



and the gauge field functions must satisfy 

w{0) = ±1 , 

for SU(2) and 

K{0) = ±2 , H{0) = , 
K{0) = , H{0) = ±2 , 
for SO (3), and the dilaton field satisfies 

0'(O) = . 



(43) 

(44) 
(45) 

(46) 



As in EYM theory |T^, [T^, it is sufficient to consider solutions with w{0) = 1 for SU(2) 
and with K{0) = 2 and H{0) = for S0(3). 
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2.5 Relations between metric and dilaton field 

Let us now derive relations between the dilaton field and the metric functions for purely 
magnetic (or electric) gauge fields of a general gauge group. For that purpose, we return 
to eq. (0) and introduce the gauge field Lagrangian La 

LM = ~\d,^d^^ + e'^''LA. (47) 

From the matter field action we then obtain the equation of motion for a static spher- 
ically symmetric dilaton field 

{ANr^^'y = -2KAe^'"^r^LA , (48) 

where for the moment we retain the dimensioned quantities. Next we return to the tt 
component of the Einstein equations, eq. (0), 

m' = 47rG QArr2$'2 - r^e^-^VL^) , (49) 

where 77 = 1 (?; = — 1) for purely magnetic (electric) gauge fields. With the help of this 
equation (|i9|), we ehminate the gauge field from the dilaton equation ( ^8[ ) 



{ANr^^y = 2r]KA\^-^-^Nr''^"'^ . (50) 

We then consider the contracted Einstein equations. For purely magnetic (or electric) 
gauge fields the contracted gauge field stress energy tensor vanishes, yielding for the 
curvature scalar 

R = SttGN^'^ . (51) 

Thus the r.h.s of the dilaton equation (|50D can be expressed purely in terms of metric 
functions 

{ANr^^J = ri^(^Am' -^Ar^R^ . (52) 

With 

Ar'^R=- [r^ {2A'N + AN')] ' + AAm' (53) 
we then obtain for the dilaton field the equation 

{ANx^'y = |7 (x^{2A'N + AN'))' , (54) 

where we have returned to dimensionless quantities. This equation can be integrated, 
yielding 

0' = |7(ln(A^iV))' + ^, (55) 



where C is an integration constant. In the following we only consider purely magnetic 
gauge fields {t] = 1). 

We are now interested in global relations between the metric and the dilaton field, 
involving the mass fi{oo) and the dilaton charge D, which we define via 



0(a;) 



D 

X 



(56) 



We further introduce the dimensionless Hawking temperature T for the metric (H) 

m B S E 

T = TsAil - 2/i')UH = TsXnAN'l^^ , (57) 

where Ts = (47ra;H)~^ is the Hawking temperature of the Schwarzschild black hole. 

We now integrate eq. (|5^) from the horizon to infinity. The asymptotic behaviour 
of the dilaton field is given by eq. (|5^). The metric functions approach constants 
asymptotically, A{oo) = 1 and A^(oo) = 1 (see sect. 2.4), while their derivatives tend 
to zero a a B 



A'{x) 



O 



1 



N'{x) 



2/i(oo) , ^ / 1 



x^ 



x-^ 



{51 



At the horizon, the metric function N{xu) vanishes, while A^'(xh) is related to the 
Hawking temperature (|57|). These boundary conditions then yield for the black holes 
the relation ||^, |30 



D = ^(fi{oo)--xlAN'\,^ 



7/i(oo) 1 - 



T 



(59) 



where /is = xi{/2 is the mass of the Schwarzschild black hole. Relation ( [59| ) holds for 
static, spherically symmetric, magnetic gauge fields of general gauge groups (provided 
the metric functions have the proper asymptotic behaviour). Thus besides the magnet- 
ically neutral SU(2) and S0(3) black holes considered here, relation ([59| ) also holds for 
the magnetically charged EMD black holes |^ and the magnetically charged SU(3) 
EYMD black holes |2|. 

The integration constant C in relation (|55|) is given by 

(60) 

Integration of 
By integrating 



C = D -7/x(oo) 



According to (|59|) it does not vanish for black holes in general ^1 
eq. (BSf) therefore does not yield a simple expression for black holes. 



Let us now consider the above relations for the regular solutions, 
eq. (|54D from zero to infinity we obtain the simple relation 



D = 7/i(oo) 



(61) 
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Consequently the integration constant C in eq. (^) vanishes for regular solutions. 
Integrating eq. (|5^) from x to infinity therefore gives the simple relation pT] 




7ln(-^«) . 



(62) 



Again, these relations are valid for static, spherically symmetric, magnetic gauge fields 
of general gauge groups. 



3 SU(2) Einstein- Yang-Mills-Dilaton Black Holes 



SU(2) EYMD theory possesses a sequence of magnetically neutral static spherically 
symmetric black hole solutions, which can be labelled by the number of nodes n of the 
gauge field function, Wn{x) ^, |^, ||, |^. These black hole solutions exist for arbitrary 
dilaton coupling 7 and for arbitrary event horizon Xh > 0. In the limit a;H 0, regular 
particle-like solutions are obtained 0, |^, |^, 0]- 

The sequences of SU(2) EYMD solutions converge to limiting solutions. In the 
following, we demonstrate numerically that in the limit n ^ 00, given a dilaton coupling 
constant 7 and an event horizon xh, the corresponding sequence of neutral SU(2) 
EYMD black hole solutions tends to the charged EMD black hole solution with the 
same dilaton coupling constant 7, the same event horizon xh and with magnetic charge 
P = 1 0. This generalizes the earlier observation 0, |], that for a given coupling 
constant 7, the sequence of regular neutral SU(2) EYMD solutions converges to the 
charged "extremal" EMD solution with the same coupling constant 7 and with magnetic 
charge P = 1. (For completeness we review the EMD black hole solutions in appendix 



Since the coupling constant 7 = 1 corresponds to the case of the low energy effective 
action of string theory, it deserves special attention. SU(2) black hole solutions for 
7 = 1 have been studied in 0, ^, where in the effect of a small dilaton mass has 
additionally been considered. 

We first discuss the properties of the black hole solutions of the stringy SU(2) 
sequence for a given event horizon. We relate the degree of convergence of the ra-th 
EYMD solution with respect to the limiting EMD solution to the location of the event 
horizon. We observe, that the smaller the event horizon xh, the larger is the value of 
the lowest n of the solution for which the functions have converged well. 

Let us begin with a small event horizon xh = 0.01. In Figs, la-e we present the 
lowest black hole solutions with an odd number of nodes, n = 1, 3, 5 and 7 for this 
event horizon. (Table 1 shows their dimensionless mass.) In particular, we present the 



A.) 



3.1 7 



= 1 
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EYMD matter functions w„ and 0„ and metric functions A^„, An and P^, where the 
charge function is obtained from the metric coefficient Qu according to 



such that P = 1 in the abehan hmiting case. For comparison we also show the corre- 
sponding functions of the hmiting EMD black hole solution, having the same horizon 
but magnetic charge P = 1. 

Fig. la shows the EYMD gauge field functions Wn and the EMD function Woo = 0. 
The value of Wn at the horizon decreases with increasing n towards the limiting value of 
zero. The location of the innermost node of Wn decreases with n, while the location of 
the outermost node of Wn increases with n, and thus Wn approaches the limiting value of 
zero in an increasingly larger region inbetween. The gauge field function wj is already 
close to zero over a region of many orders of magnitude. Obviously, with increasing n 
the EYMD gauge field functions Wn tend (nonuniformly) to the EMD function Woo = 0. 

As seen in Fig. lb, with increasing n the EYMD dilaton functions 0„ also tend to 
the EMD function (poo- The EYMD dilaton functions 0„ differ significantly from the 
EMD function (p^o only in an interior region, which decreases with n. Already for 07 
there is no notable deviation from the limiting function to be observed. 

Turning to the metric functions, we observe a similar pattern of convergence. As 
seen in Fig. Ic, the EYMD functions Nn tend to the EMD function A^oo- The limiting 
EMD function N^o rises from zero at the event horizon to a plateau A^oo ~ 1/4, and then 
rises further. Again, significant deviations between the EYMD and the EMD functions 
occur only in an interior region, which decreases with n. There the EYMD functions Nn 
develop a peak, decreasing in size with n. Again, for Nj there is no notable deviation 
from the limiting function to be observed. 

As seen in Fig. Id, the EYMD metric functions An also tend to the EMD function 
AoQ. Again, significant deviations between the EYMD and the EMD functions occur 
only in a decreasing interior region, and for Ay there is no notable deviation from the 
limiting function to be observed. 

In Fig. le the EYMD magnetic charge functions are seen and compared to 
the constant EMD magnetic charge. Poo = 1- The magnetic charge P„ tends to the 
limiting value P^o in a large inner region, roughly limited by the outermost node of 
the EYMD gauge field function Beyond this region, the magnetic charge function 
decays rapidly to zero. (Only for n = 1 there is a significant deviation from the limiting 
magnetic charge even in the inner region.) 

To understand this relation between the gauge field functions and the charge func- 
tions, we insert the asymptotic expansion (for large x) of the EYMD functions in 




(63) 
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eq. (|B3D. This yields for the magnetic charge function 

= ^K^)D'l + {wl^,^ + |/i^(oo)D2)-i + 0(1) (64) 

where the constant Wn(i) is the (l/x) expansion coefficient of the gauge field function 
Wn, Wn = (—1)" + Wn(i)/x + 0(l/x^). The value of Wn(i) is larger than /^(cxd) and D 
by six orders of magnitude, therefore the magnetic charge function has (up to a sign) 
approximately the same expansion coefficient as the gauge field function. The region 
where Wn reaches its asymptotic value is thus the same as the region where Pn decays 
to zero, i. e. the region beyond the outermost node of the gauge field function. Indeed, 
comparing the functions P^ and (w7„ + 1)^, shown in Fig. le for n = 7, we observe, that 
they are very close beyond the outermost node of the gauge field function. 

To gain more insight into the above observed pattern of convergence, let us consider 
for a moment the regular solutions. In Figs. 2a-b we show the regular EYMD gauge 
field functions Wn and metric functions for n = 1,3,5 and 7. The sequence of regular 
solutions converges to the "extremal" EMD solution which is also shown. 

Let us now discuss some features of these solutions. With increasing n, the inner- 
most node z^^-* of the EYMD gauge field functions w„ decreases exponentially to zero, as 
seen in Fig. 2a (see also Table 2), whereas the outermost node increases exponentially. 
The boundary conditions for the gauge field functions Wn do not allow for a uniform 
convergence to the EMD function Woo = 0. But the region where both EYMD and 
EMD functions are close increases exponentially with n. The EYMD metric functions 
Nn approach the value one at the origin, whereas the EMD function A^oo monotoni- 
cally approaches one quarter. Therefore there is again no uniform convergence of the 
EYMD functions to the limiting EMD function possible. But, as seen in Fig. 2b, the 
region, where the EYMD and EMD functions closely agree, again increases exponen- 
tially with n. We observe again that strong deviations of the functions from the 
limiting function occur only in an interior region, decreasing exponentially with 
n. In the vicinity of the innermost node z^^^ of w„, the function Nn has already come 
quite close to A^oo, and beyond the second node z^^-* of w„ both functions closely agree 
(with the exception of A^^i). 

The EYMD dilaton functions 0„ approach constant values close to the origin which 
decrease linearly with n. Further from the origin they approach the limiting EMD 
function 0oo- The interior regions, where the functions deviate significantly from 
the limiting solution ^oo, are limited by the respective location of the innermost node 
z^^-* of Wn (with the exception of (pi). The EYMD metric functions An exhibit an 
analogous pattern with respect to their convergence to the EMD function A^. They 
deviate significantly from the limiting function in an interior region, extending slightly 
beyond the innermost node of w„ (with the exception of Ai). 

Turning now back to the EYMD black hole solutions we observe, that the regu- 
lar solutions have considerable infiuence on the black hole solutions for small horizon 
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size and that the innermost nodes of the regular EYMD gauge field functions largely 
determine the pattern of convergence. The location of these innermost nodes z^^ is 
presented in Table 2. In particular we note, that for the black hole solutions of Figs. 1, 
the location of the innermost node of the regular gauge field function with seven nodes 
is smaller than the event horizon xu — 0.01. This fact has important consequences. 

Considering the black hole gauge field functions Wn, we see that they are large at 
the horizon as long as the horizon is small compared to the location of the innermost 
node zip of the corresponding n-th regular gauge field function. Thus Wi — are 
large at the horizon, whereas W7 is small. Consequently, wj deviates strongly from the 
limiting function only in the far region where, beyond the last node, it must approach 
its boundary value. The black hole dilaton functions 0„ deviate significantly from the 
limiting function only in the interior region between the event horizon and the location 
of the innermost node z^^^ of the corresponding n-th regular gauge field function. Thus 
there is no notable deviation left for 07 . The same is true for the metric functions A^7 
and A^. Further, we now recognize that the peak in the functions A^i — A^5, located 
shortly behind the event horizon, represents the tendency of these functions towards 
the corresponding regular metric functions. Also the functions Nn and An for n < 7 
deviate significantly from the limiting functions only in the interior region between 
the event horizon and (slightly beyond) the location of the innermost node z^^ of the 
corresponding n-th regular gauge field function. 

Choosing different values for the event horizon xh we realize that the above observa- 
tions are rather general. Considering solutions with event horizon xn = 0.1, we see that 
already the location of the innermost node of the fourth regular gauge field function is 
smaller than xn. Consequently, there are already for 04, N4 and no notable devia- 
tions from the limiting functions left. The convergence of the gauge field functions is 
somewhat slower, though is already small at the event horizon. Increasing the event 
horizon further to Xh = 1, shows that only the location of the innermost (and only) 
node of the first regular solution is still larger than the event horizon. Consequently, 
there are now already for 02, N2 and A2 no notable deviations from the limiting func- 
tions left, while ^3 is small at the event horizon. For still larger horizons, xu > 10, no 
notable deviations are left even for the functions 0i, A^i and A^. 

We demonstrate this dependence of the degree of convergence on the event horizon 
in Fig. 3, where we show the gauge field function Wj for the event horizons x-^ = 0, 0.01, 
0.1 and 1. A more quantitative analysis of the dependence of the degree of convergence 
of the functions is shown in Figs. 4a-b. There the functions Awn and A0„, defined as 
the relative deviations of the EYMD functions (l — Wn) and 0„ from the EMD functions 
(1 — Woo) and 0OO at the horizon 

AWn{Xn) = — — TT = \Wn{Xii) - Woo{Xh)\ , (65) 

[ -1- Wnn[Xv\ ) ) 
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A(/)„(xh) = r— — , (66) 

|0oo(a;H)| 

are shown as functions of the horizon for n = 1 — 5. Also shown are the locations of 
the two innermost nodes z^^ and z^^^ of the corresponding regular solutions. 

In Fig. 4a we observe, that the function Aw„(xh) crosses the line indicating the 
location of the innermost node of the n-th regular solution, zl^\ at a value of 0.27, 
for large n. Awnixn) crosses the line indicating the location of the innermost node of 
the {n — l)-th regular solution, z^li (which is located approximately halfway between 
the innermost and the second node), at a value of 0.05. (It crosses the line indicating 
the location of the innermost node of the (n + l)-th regular solution, z^li at 0.75 .) 
Considering the function A0„, shown in Fig. 4b, we see, that A^^ is much smaller 
than Awn- In particular, the function A(/)„ crosses the line indicating z^^ at a value 
of approximately 0.01 for the larger n. We observe, that for the larger n, the crossings 
lie again on straight lines, which however have a small positive slope. Concluding we 
observe, that the solutions are converged well to the limiting solution for xh S> z^\ 

Let us now consider the global properties of the solutions, the mass, the dilaton 
charge and the Hawking temperature. We recall, that these global properties are not 
independent, but related by relation (|59D. The dimensionless mass fi, the dilaton charge 
D and the Hawking temperature T/Ts for the black hole solutions with n = 1 — 7 and 
event horizons xu = 0, 0.01, 0.1 and 1 are given in Table 1. 

To discuss the inverse Hawking temperature jS = of the black hole solutions as 
a function of their mass yu(cxD), we first recall the case of the limiting EMD solution, 
which is special for 7 = 1. Here the inverse Hawking temperature is a straight line as a 
function of mass, = 47rX+ = 87ryu(oo) (see appendix A). Note, that this is the same 
relation as the one obtained for Schwarzschild black holes. However, the EMD curve 
does not extend to the origin, since the "extremal" EMD solution has a finite mass and 
a finite temperature. For P = 1 this lower limiting mass is l/-\/2- In Fig- 5 we show 
the inverse Hawking temperature as a function of the mass for the EYMD black holes 
with n = 1 — 4. As expected, the curves converge rapidly towards the limiting EMD 
curve, also shown. 

Let us now consider the convergence of the global properties of the black hole 
solutions in more detail. The mass, the dilaton charge and the Hawking temperature 
converge exponentially to the corresponding properties of the limiting EMD solution. 
To demonstrate this, we define A/i„ as the deviation of the mass of the n-th EYMD 
solution from the mass of the limiting EMD solution and make an exponential ansatz 
for the n-dependence 

A/i„ = /ioo(oo) - Hn{oo) = a^e"""" . (67) 
We define AD„ and AT„ for the dilaton charge and the Hawking temperature analo- 



15 



gously, 

AD„ = D^-Dr, = a^e-"^" , (68) 

AT„ = ^ - ^ = a^e-"-" . (69) 

The coefficients a and a still depend on the parameters, the horizon xh and the dilaton 
coupling constant 7. 

A logarithmic convergence then requires the function In (AyU„) to be a straight line 
as a function of n. This is indeed the case, as seen in Fig. 6, for the regular solutions 
and for the black hole solutions with xh > 1- However for small xh and small values 
of n, we observe deviations from straight lines. These deviations are again related to 
the relative location of the horizon to the innermost nodes z^^ of the regular solutions. 
When the location of the innermost node of the n-th regular solution becomes smaller 
than the horizon, the values of In (A/i„) for the solutions with larger n fall on a straight 
line, whereas those for smaller n are close to the straight line of the regular solutions. 

We observe, that the straight lines for the various finite values of the horizon are 
parallel, indicating the same coefficient a for all xh- (The coefficient a is different for 
the regular solutions.) For small horizons the lines are roughly equally spaced; the 
same is true for large horizons. This suggests one to parametrize the coefficients a and 
a with the following xh dependence 

«(a;H) = const. , a(xH) = 7 — rr • (70) 

Approximately, we find = 3.6, 6^ = 1.2 and 6^ = 2 for small values of the horizon, 
xr < 1, and the same but = 0.8 and 6^ = 1 for large values of the horizon, 
Xh > 1- 

The functions In (AZ)„) and In (AT„) follow a similar pattern as the function In (A/Xn) 
The coefficient a is approximately the same for all three functions, as expected, be- 
cause of relation (|59D. Also the coefficients b and 6 for the three functions are related, 
satisfying relation (^). In particular, 5^ = 5^ + 1. The constants a, b and 6 for the 
mass, the dilaton charge and the Hawking temperature obtained by a mean square fit 
to our numerical data are given in Table 3. 

In fact for large xh one obtains an analytic handle on some of the coefficients. For 
instance, it follows analytically, that 

6^ = 6d = 1, St = 2, (71) 
and, allowing for general 7, one finds 

bD = 2-fb^ , br = 2b^ , (72) 
which is indeed conffimed by the calculations. 
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3.2 7 7^ 1 

In the following we discuss the convergence of the sequences of EYMD black hole 
solutions to the corresponding limiting EMD solutions for 7 7^ 1. Note, that the values 
of 7 in 4 + ra-dimensional Kaluza-Klein theory, 7^ = (2 + n)/?T, |jl|] (with n = 1, 2 and 
4 for U(l), SU(2) and SU(3), respectively), are in the range 7 > 1. The limit 7-^0 
needs special consideration. 

The EYMD black hole solutions and their limiting EMD solutions depend smoothly 
on the dilaton coupling constant for finite 7. Not surprisingly therefore, they follow 
the same pattern of convergence for any finite 7 as for 7 = 1. However, the degree of 
convergence of an EYMD black hole solution now also depends on 7. To analyze this 
dependence, let us inspect Table 2, where the location of the innermost node 2;^^-' of the 
seven lowest regular EYMD gauge field functions is presented for the dilaton coupling 
constants 7 = 0.1, 0.5, 1 and 2. For finite 7, the sequences of nodes z^^ approach zero 
exponentially with increasing n, 

z«(7) = V^^V (73) 

The coefficients dy and Cy increase with 7. Thus the nodes zl^\'y) converge faster to 
zero for 7 > 1 and slower for 7 < 1, as compared to 7 = 1. 

The implications for the convergence of the sequences of EYMD black hole solutions 
are obvious. Let n be the number of nodes of the first EYMD black hole solution of 
a sequence with fixed horizon xh, for which the functions have largely assumed their 
limiting values. Then n becomes smaller, the larger the value of the dilaton coupling 
constant 7. This is demonstrated in Figs. 7a-b, where we show the EYMD gauge field 
function wj and the metric function A^7 of the black hole solutions with dilaton coupling 
constants 7 = 0, 0.5, 1 and 2 and event horizon xh = 0.01, as well as the corresponding 
EMD solutions P3|. According to Table 2, the location of the innermost node Zj^\'y) of 



the regular EYMD gauge field functions is larger than this horizon for 7 = and 0.5, so 
the corresponding EYMD black hole functions show large deviations from the limiting 
EMD functions in the inner regions. In contrast, for 7=1 and 2, where Zy^\'y) < xh, 
the EYMD functions have largely assumed their limiting EMD values. 

The value 7 = is special, since EYMD theory then reduces to EYM theory. The 
limiting behaviour of the sequences of EYM black hole solutions with horizon xh has 



been discussed before |T^, |T^, |I5|, ^ . For xh > 1 the limiting solutions are the RN black 
hole solutions with charge P = 1. For black hole solutions with smaller event horizons 
and for the regular solutions the limiting solutions are more complicated, having a mass 



/i(cxo) = 1, independent of the horizon ||T3|, |14|, [15], |T6|. In this case the solutions must 



be considered separately in two regions, the inner region x < 1 and the outer region 
X > 1. In the outer region the limiting solution is given by the extremal RN black hole 
solution, while in the inner region the solution is of an oscillating type with Woo 7^ 

m m m m. 
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Because of this different limiting behaviour of the black hole solutions for 7 = 
(and xh < 1) the limit 7 ^ is non-trivial. We will consider this limit separately 
elsewhere pl| . 

Let us now consider the inverse Hawking temperature (3 for general 7 beginning with 
the limiting EMD black hole solutions. For 7 > 1, the temperature in the "extremal" 
limit diverges, i. e. (3 is zero. The curve (3 versus mass then rises monotonically and 
approaches the Schwarzschild curve from below for large values of the horizon. For 
< 7 < 1 the temperature in the "extremal" limit is zero, i. e. f3 diverges in the 
"extremal" limit. The curve f5 versus mass then falls to a minimum and approaches 
the Schwarzschild curve from above for large values of the horizon. For 7 = the 
Reissner-Nordstr0m case is obtained. In contrast to the "extremal" EMD solutions, 
the extremal RN solution has a finite horizon, xn = ^ = P = 1. 

In Fig. 8, we show the convergence of the inverse temperature versus mass curves 
for the sequences of EYMD black holes with n = 1 — 4 to the corresponding curves 
for the limiting EMD black hole solutions for 7 = 0.5 and 2. (The EYM case 7 = 
is shown in Fig. 5). Depending on the dilaton coupling constant 7 and the number of 
nodes n, phase transitions can occur p, ^ |^. We observe, that for 7 > 1 no phase 



transitions occur, since the inverse temperature curves rise monotonically for all n like 
the corresponding limiting EMD curves. In contrast, for 7 < 1, phase transitions do 
occur for all members of a sequence beyond some critical Ucri'y)- The reason for this 
behaviour is, that for 7 < 1 the inverse temperature curve of the black hole solutions 
of a sequence approaches the inverse temperature curve of the corresponding limiting 
EMD black hole solutions from below. But while there is one phase transition along the 
limiting EMD curve, which starts from (3 = 00, there are two phase transitions along 
the EYMD curves (beyond some critical which start from (3 = and therefore 
develop two extrema; a maximum and a minimum. The value 7 = 1 is special in EYMD 
and EMD theory. Only for 7 < 1 phase transitions can occur. 



4 SU(3) Einstein- Yang-Mills-Dilaton Black Holes 

Let us now turn to the magnetically neutral static spherically symmetric black hole 
solutions of SU(3) EYMD theory. Besides the solutions discussed in section 3, corre- 
sponding to the SU(2) embedding, SU(3) EYMD theory possesses magnetically neutral 
static spherically symmetric black hole solutions, corresponding to the SO (3) embed- 



ding. As in EYM theory [T2|, these SO (3) solutions can be labelled by the nodes 
of both gauge field functions . Here we adopt the classification of the solutions with 
respect to the node structure of the functions (mi,M2) ^7\, |12|, 0, 

, , K{x) + H{x) , , K{x)-H{x) 

= ^ — , U2{x) = ^ ' ^ ^ . (74) 
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We denote the number of nodes of the functions Ui and U2 by rii and n2, respectively, 
and the total number of nodes by n-p. (For S0(3) solutions tit = ni + n2.) Due to 
the symmetry H{x) — >• —H{x) (see section 2.3), solutions with node structures (ni,n2) 
and {n2,ni) are equivalent, and it is sufficient to consider solutions with rii < n2- 

Previously, we have only constructed magnetically neutral SU(3) EYMD and EYM 
black hole solutions with a small number of nodes [|^ [1^. The lowest SO (3) black hole 
solution has ut = 1 and node structure (0,1), the second solution has ht = 2 and 
node structure (0,2). The third and fourth solution have also ht = 2, but both have 
the same node structure (1, 1). Therfore the node structure (ni,n2) does not seem to 
classify the S0(3) solutions uniquely. However, the third solution is a scaled SU(2) 
solution, where Ui = U2 = w, and only the fourth solution is a genuine SO (3) solution, 
where the two gauge field functions are not proportional. The two next solutions have 
riT = 3 and node structure (0,3) and (1,2), respectively. We have listed the lowest 
solutions with < 5 in Table 4, where their mass is shown for 7 = 1 and several 
values of the horizon, xn = 0, 0.5, 1, 2. 

In our previous analysis |12|, ^, we have identified two sequences of solutions: the 
genuine S0(3) sequence (0,n), containing the lowest S0(3) solution, and the sequence 
(n, n) of scaled SU(2) solutions. As in the SU(2) case, the limiting solutions of these 
two neutral sequences are magnetically charged EMD solutions, with the same event 
horizon and the same dilaton coupling constant. However, their magnetic charges are 
P = -\/3 and P = 2, respectively [0]. We discuss the convergence of these sequences in 
detail below. 

From our previous analysis 0, [T2[ it has not been clear yet, in which way S0(3) black 
hole solutions with general node structure (ni,n2) assemble to form sequences, which 
converge to limiting solutions. Below we construct sequences from the neutral SU(3) 
black hole solutions with general node structure (?T,i,n2), and we identify the limiting 
solutions of these neutral sequences with the known magnetically charged SU(3) black 
hole solutions [^, pSj] . 

4.1 (n, n) sequences 

The solutions with node structure (n, n) often come in pairs (see Table 4). One solution 
of the pair (n, n) exists always, this is the scaled SU(2) solution. The second solution 
of the pair with node structure {n,n) is a genuine S0(3) solution, whose existence 
depends on the parameters 7 and xu. 

4.1.1 Scaled SU(2) solutions 

Comparison of the equations of motion of the S0(3) embedding and those of the SU(2) 
embedding shows, that to each SU(2) solution there corresponds a "scaled SU(2)" 
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solution of the S0(3) system with precisely double the mass of its SU(2) counterpart. 
Introducing 

x = 2x , (75) 

= 2/i(x) , (76) 

K{x) = 2w{x) , H{x) = , (77) 

and 

0(x) = 0(5) , (78) 

the functions fi, w and (j) satisfy the SU(2) equations with coordinate x. Obviously, these 
scaled SU(2) solutions, having node structure (n, n), then converge to the corresponding 
EMD solutions with scaled magnetic charge P = 2. Thus these solutions and their 
properties are obtained by simply scaling the solutions of section 3. 

4.1.2 Genuine SO (3) solutions 

Although the two gauge field functions Ui and U2 of the genuine SO (3) solutions have 
the same node structure {n, n) as the scaled SU(2) solutions, they are not proportional 
to each other. We now discuss the dependence of the genuine SO (3) solutions on the 
dilaton coupling constant and on the horizon. The regular solutions exist for all values 
of the dilaton coupling constant and (supposedly) for all values of n. Considering 
the S0(3) solution with node structure (n, ra) as a function of the horizon for a fixed 
dilaton coupling constant, we observe, that it merges into the corresponding "scaled 
SU(2)" solution at a critical value of the horizon xg(7). This critical behaviour was 
first observed for vanishing dilaton coupling [|12[. In the parameter ranges where the 



genuine S0(3) solutions exist, they have a slightly higher energy than their "scaled 
SU(2)" counterparts. 

In Fig. 9, we show as an example the matter functions H and K of the genuine 
S0(3) solution (5,5) for 7 = 1 and several values of the horizon. At the critical value 
= 0.03 the solution becomes identical to the "scaled SU(2)" solution, i. e. the 
gauge field function H vanishes identically. In Fig. 10, we show the critical value of the 
horizon Xpj clS Qb function of the dilaton coupling constant 7 for the lowest odd solutions, 
= 1, 3 and 5. The critical horizon first decreases as a function of the dilaton coupling 
constant, reaches a minimum and then increases linearly with 7, = b{n)'j (with 
6(1) = 0.677, 6(3) = 0.0151, 6(5) = 0.0004) for 7 > 1. Only for 7 = the critical value 
of the horizon x^{n) converges with increasing n to a finite limit, xg(oo) = 1.3122. 

We further show in Fig. 10, the innermost nodes of the corresponding regular genuine 
S0(3) solutions as a function of the dilaton coupling constant. We observe, that for 
each value of n, the two curves closely follow each other. In general, the critical value 
of the horizon is slightly smaller than the innermost node of the corresponding regular 
S0(3) solution. Only for = 1 the two curves cross at 7 = 0.9576. (The curves 
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with the innermost nodes of the corresponding regular "scaled SU(2)" solutions also 
have the same shape but lie somewhat higher.) From our previous considerations on 
SU(2) black holes, this behaviour is not unexpected. There, when the horizon of the 
black hole becomes larger than the innermost node of the regular solution, the solution 
loses most of its structure. Here, since the genuine SO (3) black hole solutions have 
more structure than the "scaled SU(2)" black hole solutions, this additional structure 
disappears, and thus the genuine SO (3) solution disappears, when the horizon comes 
close to the vicinity of the innermost node of the regular SO (3) solution. 

This observation demonstrates even more the importance of the innermost nodes of 
the regular solutions for the black hole solutions. For the SU(2) solutions, the innermost 
nodes of the regular solutions rule the pattern of convergence of the sequences of black 
hole solutions. For the genuine (n, n) S0(3) solutions, they determine the range of 
existence of the corresponding black hole solutions. 

4.2 (0, n) sequences 

The two lowest neutral black hole solutions of the SO (3) embedding have node structure 
(0, 1) and (0,2). They are the first members of the sequence (0,?t,). For a given event 
horizon and dilaton coupling constant, this sequence tends to the EMD black hole 
solution with the same event horizon, the same dilaton coupling constant and with 
magnetic charge P = \^ . Table 5 shows the dimensionless mass, the dilaton charge 
and the Hawking temperature for the lowest solutions of this sequence for 7 = 1 and 
horizon xh = 1. 

Previously, we have only demonstrated the convergence of the global properties of 
this sequence. Here we consider in addition, the convergence of the spatial functions 
of the solutions of this sequence with increasing node number n. As in the SU(2) case, 
the location of the innermost node of the regular solution with n nodes determines the 
degree of convergence of the corresponding black hole solution with n nodes and horizon 
xh- (Remember, that the nodes of the function U2{x) correspond to the intersections 
of the functions K{x) and H{x).) We observe that the black hole solution with n 
nodes has already largely converged to the corresponding limiting EMD solution, if the 
innermost node of the regular solution with n nodes is smaller than the horizon. It has 
not yet converged well, if the innermost node of the regular solution is larger than the 
horizon. Thus the pattern of convergence of the (0, n) sequence is fully analogous to 
the SU(2) case. We demonstrate this in the following with a few examples. 

Let us begin again with the case 7 = 1, the case of the low energy effective action of 
string theory. In Figs, lla-c, we show the matter functions Kn and Hn, and the metric 
functions and for the small horizon xh = 0.2 for n = 1, 3 and 5. (The matter 
functions 0„ and the metric functions are analogous to the SU(2) case and contain 
no relevant new information.) From Fig. 11a, we observe, that with increasing n, both 
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EYMD gauge field functions Kn and Hn approacli tlie same limiting function 

i^n(a:)™^, (79) 

where / is constant 

f{x) ^ 1 . (80) 

To relate this limiting function / to the corresponding EMD solution with charge 
P = a/S we inspect the EYMD equations of motion for the case K = H = f/\/2. 
(These equations are presented in appendix B.) Indeed, for / = 1 these equations 
reduce to the EMD equations of motion for a black hole with charge P = ^/S. 

Let us now consider the degree of convergence of the black hole solutions of Figs. 11a- 
c with xh = 0.2. For that purpose, we inspect the location of the innermost nodes of 
the gauge field functions of the corresponding regular solutions, which are given in 
Table 6. As for the regular SU(2) solutions, with increasing n the innermost nodes 
decrease exponentially to zero. For the example of Figs, lla-c, the innermost nodes of 
the regular solutions are larger than the horizon for n = 1 and 3 and smaller than the 
horizon for n = 5. We conclude, that the solution with n = 5 should be close to the 
limiting EMD solution, also shown in Figs, lla-c, whereas the solutions with n = 3 
and n = 1 should still show large deviations from the limiting solution. This is indeed 
the case. The functions 05, and are almost indiscernible from the limiting EMD 
functions, whereas the functions 03, and ^3 deviate significantly from the limiting 
functions in the region interior to the innermost node. The gauge field functions Kn 
and Hn converge more slowly (and non-uniformly because of the boundary conditions 
at infinity), but ^ i^s ^ l/\/2 already in a large region of space. Also the charge 
function P^ has assumed its limiting value = 3 in a large region. In contrast, the 
charge function P^ still deviates slightly from the limiting value in the inner region and 
Pi deviates strongly. 

We demonstrate the dependence of the black hole solutions on the horizon xh in 
Fig. 12. There the gauge field functions and Hr, are shown for the event horizons 
xh = 0, 0.5, 1 and 2. For the larger horizons, even the gauge field functions have 
converged well in the inner region. (For the regular solutions the boundary conditions 
inhibit a uniform convergence at the origin.) As for the SU(2) solutions, a more quan- 
titative analysis of the degree of convergence of the solutions is obtained by studying 
the functions AKn, AHn and A(f)n, where the functions AKn and AHn, are defined 
as relative differences analogously to A0„ (eq. (^)). We observe, that the function 
AKn{xi{) crosses the line indicating the location of the innermost node of the n-th 
regular solution, zl^\ at the value 0.75 for the larger n, while the function AHnixn) 
crosses at the value 0.60. Both functions have fallen to the value 0.05, when they cross 
the line indicating the location of the second node of the n-th regular solution, z^^ . The 
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function A0„ is again much smaller than Aii'„ or AHn and crosses the line indicating 
zl^^ at a value on the order of 0.05 for the larger n. 

The SO (3) solutions depend smoothly on the dilaton coupling constant 7. The 7 
dependence of the convergence properties of the sequences of black hole solutions is 
analogous to the SU(2) case. With increasing 7, the convergence of the sequences is 
faster, with decreasing 7 the convergence is slower. This is seen from Table 6, where 
the innermost nodes of the lowest regular solutions of the (0, n) sequence are given for 
7 = 0.1, 0.5, 1 and 2. 

Again, the limit 7 — needs special consideration and will be discussed elsewhere 
33. However, we already note, that for 7 = the limiting solution is more complicated 



and analogous to the limiting SU(2) solution. For xu > P (P = V^), the limiting 
solution is the RN solution, while for xu < P the limiting solution is the extremal RN 
solution for x > P and apparently an oscillating solution for x < P. 

Let us now consider the convergence of the global properties of the solutions. Turn- 
ing to the thermal properties of the black hole solutions first, we show in Fig. 13 the 
inverse Hawking temperature /3 as a function of the mass for the EYMD black holes 
for n = 1 — 4 and 7 = 1 as well as 7 = 0. As in the SU(2) case, the EYMD inverse 
temperature curves converge rapidly towards the limiting EMD and RN inverse tem- 
perature curves. Consequently, we also observe phase transitions, depending on the 
dilaton coupling constant and on the number of nodes. For instance, for 7 = and 
n = 1, there are two critical values of the mass, fii = 1.561 and fi2 = 1.588. The critical 
behaviour disappears beyond 7 = 0.0293 for the lowest solution, while it occurs up to 
larger values of 7 for the solutions with larger n. For 7 > 1 no phase transitions occur. 

The convergence of mass, dilaton charge and Hawking temperature is exponential, 
analogous to the SU(2) case. For the dilaton coupling constant 7 = 1, these global 
properties are given in Table 5 for n = 1 — 5 together with their limiting values and the 



constants of the exponential approximation formulae eqs. (|67D-(|69D [p6l . For large val- 
ues of the horizon, the Xn-dependence in the approximation formulae can be extracted 
according to eq. (|70D . The corresponding parameters are shown in Table 7. 



4.3 (1, 1 + n) sequences 

Let us now address the question, in which way the solutions with general node structure 
(ni,n2) assemble into sequences and to what limiting solutions these sequences con- 
verge. The existence of the sequence (0, n), where one gauge field function has no node, 
whereas the other function has an increasing number of nodes, suggests one to examine 
the set of solutions, where one of the gauge field functions has one node, whereas the 
other function has an increasing number of nodes. We thus consider the set of solutions 
(1, 1 + n). The first solution of this set is the solution (1, 2); it has n = 1. (The solu- 
tion (1,0) is equivalent to the solution (0, 1) and belongs to the (0,n) sequence. The 
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solutions (1, 1) are members of the (n, n) sequences.) In the following, we demonstrate 
that for any value of the dilaton coupling constant 7 and of the horizon Xh, the set of 
solutions (1, 1 + n) indeed forms a sequence, tending to a limiting solution. 

To identify the limiting solution, let us consider an example. In Figs. 14a-c we 
show the matter functions Kn and and the metric functions Nn and for the 
horizon xh = 0.2 and = 1, 3 and 5 for 7 = 1. As for the (0, n) sequence, in the limit 
n ^ 00 both EYMD gauge field functions and approach a single function /(x). 
However, in contrast to the (0, n) sequence, this function is not a constant; it depends on 
X and has one node. Obviously the limiting solution cannot be a magnetically charged 
EMD black hole. Instead the limiting solution turns out to be the lowest magnetically 



charged SU(3) black hole solution found in |26]. This limiting solution is also shown in 
Figs. 14a-c. 

In fact, the simpler system of equations obtained with only one gauge field function, 
(see appendix B) 

K{x) = H{x) = f{x)/V2 , (81) 
possesses SU(3) black hole solutions with magnetic charge P = ^/S. Such charged black 



hole solutions have been constructed before for the dilaton coupling constants 7 = 1 p6 
and 7 = |^|. The equations for these charged black hole solutions closely resemble 
the SU(2) equations (and correspond to the SU(2) x U(l) equations). Not surprisingly 
therefore, the charged black hole solutions form a sequence, which can be labelled by 
the number of nodes j of the single gauge field function / p6|, Like the EMD 



black holes, for finite dilaton coupling constant 7, these charged black hole solutions 
with j nodes exist down to arbitrary small values of the horizon. In the "extremal" 
limit xh = 0, the solutions present naked singularties like their EMD counterparts. 
(The "extremal" limit is discussed in appendix B.) In contrast, for vanishing coupling 
constant 7, the extremal black hole solutions have a finite horizon, xn = P (P = a/3), 
analogous to the RN case [p5|] . 

Turning back to the (1, 1 +n) sequence of EYMD solutions, we observe, that for any 
finite value of the dilaton coupling constant and any finite event horizon, the sequence 
of neutral black hole solutions tends to the magnetically charged SU(3) solution with 
one node, j = 1. This limiting solution has the same dilaton coupling constant, the 
same horizon, and it has magnetic charge P = a/3. For finite dilaton coupling constant 
7, the "extremal" solution with j = 1 represents the limiting solution of the regular 
neutral S0(3) sequence (1, 1 + n). However, for zero coupling constant 7, the extremal 
black hole solution has a finite horizon, xh = P, therefore the limiting solution is again 
more complicated for values of the horizon smaller than P . 



The convergence properties of the solutions of the (1, 1 + n) sequence are again 
analogous to those of the S0(3) {0,n) sequence and the SU(2) sequence. The location 
of the horizon with respect to the innermost nodes of the function U2 of the regular 
solutions determines the degree of convergence of the black hole solutions. The inner- 
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most nodes of the regular solutions are shown in Table 6 for several values of the dilaton 
coupling constant 7. For the example of Figs. 14a-c, the horizon is slightly smaller than 
the innermost node of the regular solution with n = 3. Consequently, the black hole 
solution with n = 5 has already converged well, the solution with n = 3 has almost 
converged, and the solution with = 1 is far from convergence. Note the interesting 
shape of the limiting charge function P^{x), which is already approached well in the 
inner region by P5 and in a smaller region by P^. The limiting charge function has two 
plateaus, an inner plateau at 3.790 and an outer plateau at the magnetic charge of the 
solution, = 3. 

The global properties of the lowest solutions of the (1, 1 + n) sequence for 7 = 1 and 
xh = 1 are given in Table 5. The inverse temperature as a function of mass is shown 
for the lowest solutions of the (1, 1 + n) sequence and for the charged limiting j = I 
solution in Fig. 15 for 7 = 1 and 7 = 0. The exponential convergence of the global 
properties is demonstrated in Table 5 and Table 7 for 7 = 1. 



4.4 {j,j-\-n) sequences 

From the above results, it is clear, in which way the remaining SO (3) solutions fall 
into sequences. The next sequence is formed by the solutions (2, 2 + n), beginning with 
the solution (2,3). In Figs. 16a-b, we show the matter functions Kn and and the 
charge functions for n = 1, 3 and 5, the horizon xh = 0.2, and the dilaton coupling 
constant 7 = 1. In the limit n ^ 00, both EYMD gauge field functions Kn and if„ 
approach a single function f{x) which has two nodes. This limiting function is the 
gauge field function of the magnetically charged SU(3) EYMD black hole solution with 
two nodes, j = 2, and charge P = a/3. 

In general, we observe, that for any finite dilaton coupling constant 7, the (2, 2 + n) 
sequence of EYMD solutions tends to the corresponding magnetically charged SU(3) 
solution with two nodes, j = 2, and magnetic charge P = ^/3. Considering the conver- 
gence properties of the solutions of the (2, 2 + n) sequence, we observe that they are 
again analogous to those of the above sequences (0,n) and (1, 1 + n). The location of 
the horizon with respect to the innermost nodes of the regular solutions determines the 
degree of convergence of the black hole solutions. 

In the example of Figs. 16a-b, the horizon is located well behind the innermost node 
of the regular solution with n = 3, as can be seen in Table 6 Consequently, the 



black hole solutions with n = 5 and n = 3 have already converged well, whereas the 
solution with n = 1 is still far from convergence. In particular, also the charge function 
P^ of the limiting solution is well approached in the inner region by P^ and P^. We 
note, that the inner plateau of the limiting charge function has increased in height to 
3.964 and has also increased in size. We further note, that the deviations AKn fall on 
top of each other for j = 0, 1 and 2, when regarded as functions of the scaled horizon 



25 



xh/^^s (i)' the same is true for the deviations AiJ„. 

The global properties of the lowest solutions of the (2, 2 + n) sequence are given in 
Table 5. The exponential convergence with n of the global properties is demonstrated 
in Tables 5 and 7, where we anticipate also an exponential convergence with j. 

By generalizing the above observations we see, that (apart from the (n, n) solutions 
considered separately above) the neutral SO (3) EYMD black hole solutions assemble 
into sequences {j,j + n), with fixed index j > and running index n > 1. For any 
finite value of the dilaton coupling constant and any finite event horizon, the sequence 
(j, j+n) of neutral black hole solutions tends to the magnetically charged SU(3) solution 
with j nodes, the same dilaton coupling constant, the same horizon, and with magnetic 
charge P = ^/S. For finite dilaton coupling constant 7, the "extremal" solution with j 
nodes represents the limiting solution of the regular neutral S0(3) sequence (j, j + n). 
(For coupling constant 7 = 0, the limiting solution is the RN solution if xh > P] it is 
again more complicated, if xh < -P pi|.) 

Let us finally consider the sequence of the sequences (j, j + n) of SO (3) solutions, 
with running index j. The lowest SO (3) sequence has j = 0. It approaches an EMD 
solution with magnetic charge P = -\/3 since / = 1. The (j, j + n) sequences approach 
the EYMD solutions with magnetic charge P = a/3 and j nodes. In the limit j ^ 00, 
the limiting solutions of the sequences tend to a limiting solution themselves. This 
limiting solution has / = 0, and therefore it is again an EMD solution. But this 
limiting EMD solution has magnetic charge P = 2. (It is also identical to the limiting 
solution of the sequence of scaled SU(2) solutions.) The inverse temperature curves of 
the limiting charged solutions for j = 1 — 4 are shown in Fig. 17, together with the 
curves for j = and j = 00 for 7 = 1 and 7 = 0. 

5 Conclusion 

We have considered sequences of magnetically neutral, static, spherically symmetric 
black hole solutions of SU(3) Einstein- Yang-Mills-dilaton theory, based on the SU(2) 
embedding and on the SO (3) embedding. Such sequences of black hole solutions exist 
for arbitrary dilaton coupling constant 7 and arbitrary event horizon Xh. In the limit of 
vanishing event horizon, sequences of regular EYMD solutions are found. In the limit 
of vanishing dilaton coupling constant, sequences of EYM solutions are obtained. 

The equations of motion for the EYMD system, together with the boundary condi- 
tions, allow for relations between the metric and the dilaton field. Black hole solutions 
satisfy a relation between mass and dilaton charge, D = '-f{fi{oo) — 27ix^T), where T is 
the Hawking temperature. Regular solutions satisfy two simple relations, D = 7/^(0x3) 
and (f){x) = 7 ln{^y—gtt)■ These relations hold for static, spherically symmetric solutions 
with magnetic gauge fields of general gauge groups. 

The members of the SU(2) EYMD sequences are labelled by the number of nodes 
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n of the single gauge field function. In the limit of large n, for a given dilaton coupling 
constant 7 and a given horizon xh, the corresponding sequence of black hole solutions 
tends to the EMD black hole solution with magnetic charge P = 1, the same 
coupling constant and the same horizon. The corresponding sequence of regular EYMD 
solutions tends to the corresponding "extremal" EMD solution. 

The solutions of the SO (3) embedding can be labelled by their node structure 
(^1, 712), representing the number of nodes of both gauge field functions. The S0(3) so- 
lutions fall into sequences of two types. Most solutions belong to the first type, having 
node structure (j, j + n), where j labels the sequence itself and n labels the members 
of the sequence, i. e. the lowest sequence has j = 0, the next j = 1, etc. The solutions 
of the second type have node structure (n, n). Beside the "scaled SU(2)" solutions, 
there exist also genuine SO (3) solutions of this type. With regard to the genuine SO (3) 
solutions, the above labelling of the solutions by their node structure appears to be 
unique. 

In the limit of large n, for a given dilaton coupling constant 7 and a given horizon 
xu, the corresponding sequence of black hole solutions with node structure (j, j + n) 
and finite j tends to the magnetically charged SU(3) EYMD black hole solution with j 
nodes pGf, which has the same coupling constant, the same horizon and magnetic charge 
P = a/3. The correspoding sequences of regular solutions tend to the corresponding 
"extremal" magnetically charged EYMD solutions. 

For j = and j — * 00, the magnetically charged EYMD SU(3) black hole solutions 
have constant gauge field functions. Therefore the limiting solutions in these cases are 
again charged EMD black hole solutions, having magnetic charges P = ^/3 and P = 2, 
respectively. 

The genuine S0(3) solutions of type {n,n) differ from those of type + n) in 
an important aspect. They do not exist for all values of the dilaton coupling constant 
7 and the horizon xh- Instead, for each value of the coupling constant, there exists a 
critical value of the horizon, where the genuine S0(3) solution {n,n) merges into the 
corresponding "scaled SU(2)" solution. 

In the case of vanishing dilaton coupling constant, the sequences of neutral EYM 
black hole solutions tend with increasing n to magnetically charged EYM black hole 
solutions and to RN black hole solutions (for (0,n) and scaled SU(2) {n,n)), as 
long as the horizon is larger than the magnetic charge. If the horizon is smaller than 
the magnetic charge, the limiting solutions are more complicated [|T^, |T^, |T^, 

All sequences of black hole solutions exhibit the same pattern of convergence. The 
degree of convergence of a given EYMD black hole solution, labelled by n within the 
corresponding sequence, depends on the relative location of the horizon to the innermost 
node of the corresponding regular EYMD solution. If the horizon is larger than the 
location of this node, then the functions of the EYMD black hole solution have largely 
converged to the functions of the limiting solution. If the horizon is smaller than the 
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location of this innermost node, then the functions of the EYMD black hole solution 
still differ appreciably from those of the limiting solution, with a tendency to approach 
the functions of the regular EYMD solutions for small x. 

The location of the innermost nodes of the regular EYMD solutions converges to 
zero exponentially, for all sequences. The larger (smaller) the dilaton coupling constant, 
the faster (slower) the convergence. For 7 = the location of the innermost nodes 
converges to a finite value for the various sequences. 

For all sequences, the convergence of the global properties of the EYMD solutions, 
such as mass, dilaton charge and Hawking temperature, to those of the limiting charged 
solutions is exponential. The coefficients depend on the dilaton coupling constant 7, 
yielding a fast convergence for large 7 and a slow convergence for small 7. 

Obviously, the thermodynamic properties of the EYMD black hole solutions also 
tend to those of the corresponding limiting solutions for large n. The occurrence of 
phase transitions, where the specific heat changes sign, depends on the dilaton coupling 
constant 7. Only for 7 < 1 phase transitions can occur. 

Addressing the stability of the solutions, we note that the SU(2) and SU(3) EYM 
black hole solutions are unstable [0, |T|, |21|, g |2|, g, and so are the SU(2) EYMD 



black hole solutions [Q. We therefore conjecture, that the genuine S0(3) EYMD black 
hole solutions are also unstable. It appears to be interesting to study the number of 
unstable modes of the genuine SO (3) black hole solutions and look for a relation to the 
number of nodes of the solutions, in analogy to the SU(2) case, where the black hole 
solution with n nodes has 2n unstable modes [^, [21[| . 
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x = x[i-:^Y^' , (82) 



6 Appendix A: Einstein- Maxwell-Dilaton Black Holes 

We here briefiy recall the well-known static, spherically symmetric EMD black hole 
solutions with magnetic charge P Following the notation of we introduce the 

coordinate X via 

X 
~X 

X = X_ then corresponds to the origin x = 0. (The coordinates x and X correspond 
to R and r of 0, respectively.) The metric then takes the form 

ds^ = -X^dt^ + X~^dX^ + x^de^ + sin^ ed(j)^) , (83) 

with 



X / V X 
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The non-vanishing component of the Maxwell field strength tensor is given by 

Fe^ = Psme (85) 

and the dilaton field by 

, Y .27/(1+72) 
e2'/'=(l_4_) . (86) 



X 

The parameters X+ and X_ are determined by the regular event horizon xh 



XH = X+ 1-— (87) 



and by the magnetic charge P 

'X.XJ 



P 



V 1 + 7 V 

The black hole solutions have mass /^(oo) 

1 + 7 

and dilaton charge D 



D = ■ (90) 



7X_ 

Thus they satisfy a quadratic relation between mass, dilaton charge and magnetic 
charge 

fi'ioo) + D'-P'=( • (91) 



2 

The EMD black hole solutions exist for arbitrary event horizon xb_. The limit 
xh — > is known as "extremal" limit. The "extremal" solutions have naked singularities 
at the origin. For finite horizon the solutions satisfy relation (^9]), whereas in the 
"extremal" limit they satisfy relations (|6lD and (|62l). For "extremal" EMD solutions 
the quadratic relation (pT| ) between mass, dilaton charge and magnetic charge simplifies 
to 



/x = P/Vl + 7'. (92) 
For comparison let us also note the metric functions A and N of the metric (P) 

For the "extremal" solutions 

^(0) = (ItW • (94) 
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6.1 7 = 1 

For 7 = 1 the above formulae simplify considerably. The relation between the coordi- 
nates X and X is easily inverted 



X_ + J XI + 4x2 
X = . (95) 



Eliminating the parameter X_ via 



^- = ^ , (96) 



the parameter X+ is determined by 

X+ = ^xl + 2P2 (97) 

for a given horizon xh and magnetic charge P. The mass, dilaton charge and Hawking 
temperature then are 

/x(oo) = X+/2 , (98) 

D = PVX+ = , (99) 

//i(ooj 



and 



6.2 7 = 



T/Ts = xh/X^ = . (100) 

//i(OOj 



For 7 = the dilaton field decouples. The black hole solutions with magnetic charge 
P are simply the Reissner-Nordstr0m solutions with metric (H), where 



X=(l-=^ + — ), A = l . (101) 



2u P 
+ - 

The non- vanishing component of the Maxwell field strength tensor is given by eq. (|85| ) as 
above. For magnetically neutral black holes the metric (|101|) reduces to the Schwarzschild 
metric, where P = 0. 



7 Appendix B: Charged SU(3) Einstein- Yang-Mills 
Black Holes 

A restricted subset of SU(3) EYMD black hole solutions with K = H has been studied 
previously |^ . These solutions correspond to magnetically charged black holes . 
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To derive their equations of motion, let us put K = H = m. the SO (3) equations, 
eqs. (P^)-(P(]|). This yields for the metric functions 



. \2 3 



X 



2 ^'2 ^ g270 fl2 



A . 



For the matter field functions we obtain the equations 

{e^^'t'ANfy = e^^^^Af [f - l) 



{ANx'(j)y = 2'yAe 



2x^ 



3 

2^2 



(102) 
(103) 

(104) 
(105) 



These equations differ from the SU(2) equations only by the presence of the term ^ in 
the equations of /i and 0. In fact they correspond to magnetic SU(2) x U(l) equations. 
In general a term of the form ^ is present for magnetically charged solutions with 



charge P. For the asymptotic behaviour of the gauge field function 

fix) 1 , 

the above solutions have magnetic charge P = y^, while for 

fix) ™ , 



(106) 



(107) 



their magnetic charge is P = 2. 

Obviously, the condition i^' = if is in conflict with the boundary conditions (^) 
and (|3^), imposed on the gauge field functions of neutral black hole solutions at infinity. 
Neither can the boundary conditions (|4^) and (^Sf) , imposed for the gauge field function 
of regular solutions at the origin, be satisfied for solutions with K = H. Therefore these 
solutions have no regular limit . For 7 = there exists for any n an extremal black 
hole solution with finite horizon xh = P, in close analogy to the extremal RN solution 
25| . For 7 7^ the black hole solutions exist for any n down to arbitrary small horizon. 
In the limit Xh — * an "extremal" solution is obtained ^ . 

In the following we show, that this "extremal" solution is closely analogous to the 
"extremal" EMD solution. To study the limit a; —> 0, let us introduce new functions G 
and r by 

xG' 



G 



7X0' 



G 



with 



G^Gox , r ^ 1 



108 



(109) 
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Inserting the small x expansion for the functions G and / in eqs. ( p.02| ), ( |103| ) and ( |105| ), 
we obtain to lowest order 

f' = ~ + \<^,P' . (110) 

(111) 

and 

(ANx)' = -f^AGlP^ . (112) 

There are no contributions from the gauge field function to this order. Solving eqs. ( |11(J| ) 
and ( |1 1 1| ) gives to lowest order in x 

^=lTTV'^V'+^on-, (113) 

A = Aqx^ , (114) 
where Aq is an integration constant. Eq. ( |112| ) then leads to 

GlP' = ^ . (115) 



Thus we find to lowest order in x 



1 l + V 



G= . ^ =3: , (117) 



and at the origin 

iv(o) = i-^r . (118) 

To determine Aq we note, that eq. ( |0^ holds also for the "extremal" EYMD solutions. 
Rewriting eq. (|62[) in terms of G(x) yields 

G = {A^Ny'/'' . (119) 

Inserting the small x expansions, eqs. ( |1 14|) , (|117|) and ( [118| ), we find the constant Aq, 

Aq = 1^ {^i + j2py^ . (120) 

Go, ^0 and iV(0) coincide with the corresponding expressions for the "extremal" EMD 
solutions in the limit X. = X_. 
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Using the functions G and F instead of and 0', we have constructed the "extremal" 
solutions numerically with the boundary conditions 



G™ 1 , r(0) = 1 , (121) 

and the boundary conditions of the regular solutions for the other functions. 

Let us finally calculate the inverse Hawking temperature of the "extremal" solutions. 
At the horizon the equations for ^ and F are given by 

^^n = Gl[ ]h, rH = ^^[ ]h. (122) 

(The index H indicates the value of the functions at the horizon, [ ]h abbreviates an 
irrelevant expression.) Eliminating GW ]h in eq. (|122|) we find 



Mi = ^if^* , (123) 



In the limit — with Fh — >■ 1 and /ih = ^ we obtain 

(Note, that the limit x — for the "extremal" solutions coincides with the limit xh — ^ 
of the black hole solutions for the functions G, F, A, /, but not for A^.) 

Using the small x behaviour of A, eqs. (|114|) and (|12CI|) , the inverse Hawking tem- 
perature becomes in the limit of vanishing horizon 

1 .,2 1 



/3(xh) 4vr l^^l+^^py . (125) 

Thus for the "extremal" EYMD solutions (3 is finite for 7 = 1, it diverges for 7 < 1 
and vanishes for 7 > 1, completely analogous to the inverse Hawking temperature of 
the "extremal" EMD solutions. 

The limit of small horizon for the EYMD black hole solutions on the one hand and 
the limit of small x for the "extremal" EYMD solutions on the other hand coincide 
with the limiting behaviour of the corresponding EMD solutions. The reason is clearly, 
that there is no explicit contribution from the gauge field in this limit. Therefore, for 
7 = 1, /3 has the same value for the "extremal" EMD solution and for the "extremal" 
EYMD solutions with j nodes, as seen in Fig. 17. 
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[31] An exception are the "extremal" magnetically charged EMD and EYMD solutions 
(see Appendices A and B). 

[32] The relations ( |6TD and ( |62D generalize the relations obtained previously for SU(2) 
and 7 = 1 [0, ||, ||, |2^. Ref. contains a scaled version of relation (^) for SU(2) 
and general 7. 

[33] Relation (^) for can be generalized for 7 7^ 1. In order to define the charge 
function for 7 7^ 0, we start from —x^g^^^"^ = {X — X_){X — X^). We introduce 
the variable ^(X) = X--f^^ = X--fD, and employ + = 2(/i(cx>) +7L') 
and X_ = ^D. We then find -x^g^f^^^ = + ^2 _ ^2 _ 2^/1(00), and we 
define the charge function P{x) for the EYMD system by 

P\x) = -x^gu - eix) + 2e(x)/i(oo) + , 

with the function ^(x) implicitly given. With y"^ = — the equation for the 
charge function can be put into the form —x'^gu = P^{y) — 2/i(cx3)v^?/^ + D'^ , 
which coincides with ( |63D for 7 = 1 (since y"^ = x^). 

[34] B. Kleihaus, J. Kunz and A. Sood, in preparation. 

[35] The coupling of ref. ^ is related to ours via k = V2'j. 

[36] The masses /i(oo) of the regular solutions are well approximated by the formula 

/i(oo,7,n) = [v^- (7r/2)exp(-4n/3)]/(v/rT^) 

[37] For 7 = 1 and xh < 1 many data points of the sequences (0,n), (1, 1 + n) and 
(2, 2 + n) for n < 5 do not yet fall on straight parallel lines. This makes it difficult 
to extract the asymptotic behaviour for large n. Employing 

A/i„,,(oo) =/ioo,,(oo) -^„,,(oo) = %e-("''"+^-^) , 

(analogously for D) and 

AT,,, = (T/T,)„,, - (T/T,)„o,, = :^e-(°-'^+^-^') , 

we find for the coefficients ~ 6.7, bo ~ 19.3 and 6t ~ 25.6, and for the 
exponents ~ an ~ ay ~ 2.8 and ~ q;d ~ ar ~ 3.7. 

[38] Table 6 shows, that the innermost nodes of the regular solutions decrease approx- 
imately exponentially with j for fixed n. 

[39] The equations of [^, ^ contain a scale factor Rg. They agree with our equations 
(|102|) - (|105|) for Re = 2. The black hole solutions of [|^ are presented for R^ = 1, 



and therefore tend to a scaled limiting solution with charge P = 1. 
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6 




.70709 


.70708 


.70535 


.57735 


7 




.70710 


.70709 


.70535 


.57735 


T/T5 


n 




.0 


.01 


.1 


1.0 


1 




.39365 


.39601 


.41750 


.64915 


2 




.16649 


.16900 


.19299 


.58048 


3 




.06773 


.07028 


.09870 


.57744 


4 




.02738 


.03000 


.07256 


.57735 


5 




.01106 


.01388 


.07059 


.57735 


6 




.00446 


.00809 


.07054 


.57735 


7 




.00180 


.00711 


.07053 


.57735 



Table 1 

The mass /i(oo), dilaton charge D and Hawking temperature T/Ts for the SU(2) EYMD 
regular solutions (xh = 0) and black hole solutions (xh = .01, .1 and 1) with up to 
seven nodes n and dilaton coupling constant 7 = 1. 
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n 


7 




0.1 


0.5 


1.0 


2.0 


1 


1.53846 


1.42382 


1.35402 


1.62119 


2 


1.08289 


.80649 


.53637 


.44233 


3 


.95933 


.52922 


.21751 


.10992 


4 


.92289 


.36395 


.08791 


.02606 


5 


.90328 


.25272 


.03550 


.00612 


6 


.88668 


.17577 


.01434 


.00143 


7 


.87082 


.12229 


.00579 


.00034 




.018 


.363 


.907 


1.451 


dj 


.988 


1.550 


3.308 


8.661 



Table 2 

The location of the innermost node -2^^-* (7) for the SU(2) EYMD regular solutions up 
to seven nodes n for several values of the dilaton coupling constant 7. The nodes are 
well approximated by z^^^'y) = d^e'^''", with constants d^, also shown. 





< xh < 1 
/i(oo) D T/Ts 


10 < xh < cxo 
^(00) D T/Ts 


s 
h 

a 


1.916 1.951 2.941 
1.156 2.960 3.361 

3.623 3.591 3.568 


1.001 1.002 2.006 
0.782 1.594 1.607 

3.591 3.630 3.624 



Table 3 

The constants 5, 6, a of the least square fits for A//„ = //oo(oo) — iin{oo) — , \sa ^~°''^^ 
(and analogous for D and T /Tg ) for the SU(2) EYMD sequence for 7 = 1. The fits 
have been done separately for the regions < xh < 1 and 10 < xh < 00. 
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UJ + n) 








0. 




0.5 




1.0 




2.0 


(0,1) 


.90853 


1 


02708 


1 


16497 


1 


49442 


(0,2) 


1 


14153 


1 


20944 


1 


30297 


1 


57396 


(1,1)* 


1 


15397 


1 


25989 


1 


38135 


1 


67341 


(1,1) 


1 


18797 


1 


27399 


1 


38437 






(0,3) 


1 


20383 


1 


24529 


1 


32140 


1 


58069 


(0,4) 


1 


21958 


1 


24966 


1 


32279 


1 


58111 


(0,5) 


1 


22347 


1 


24998 


1 


32287 


1 


58114 


(1,2) 


1 


32130 


1 


37840 


1 


46227 


1 


71399 


(1,3) 


1 


36195 


1 


39923 


1 


47209 


1 


71729 


(2,2)* 


1 


36967 


1 


42005 


1 


49382 


1 


73014 


(1,4) 


1 


37300 


1 


40151 


1 


47278 


1 


71750 


(2,2) 


1 


37417 














(2,3) 


1 


39866 


1 


43109 


1 


49813 


1 


73148 



Table 4 

The mass of the lowest regular (xh = 0) and black hole (xh = 0.5, 1 and 2) S0(3) 
EYMD solutions with a total number of nodes < 5 and dilaton coupling constant 
7 = 1. The asterix indicates the scaled SU(2) solutions. Note that the genuine S0(3) 
black hole solutions with node structure (n, n) do not exist for all horizons. 
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UJ + n) 

(0,1) 
(0,2) 
(0,3) 
(0,4) 
(0,5) 
(0,oo) 



1.16497 
1.30297 
1.32140 
1.32279 
1.32287 
1.32288 



D 
.87027 
1.08770 
1.12979 
1.13364 
1.13388 
1.13389 



T/Ts 

.58942 

.43054 

.38322 

.37830 

.37798 

.37796 



a 



(1,2) 
(1,3) 
(1,4) 

(1,5) 
(1,6) 
fl,oo) 



2.798 
6.529 



2.786 
17.582 



1.46227 
1.47209 
1.47278 
1.47283 
1.47283 
1.47283 



1.26493 
1.28839 
1.29029 

1.29041 
1.29042 
1.29042 



2.780 
22.235 



.39467 
.36740 
.36499 

.36484 
.36483 
.36483 



a 



(2,3) 

(2,4) 
(2,5) 
(2,6) 
(2,7) 
(2,oo) 



2.773 
.180 



2.748 
.443 



1.49813 

1.49843 
1.49845 
1.49845 
1.49845 
1.49845 



1.32829 

1.32917 
1.32922 
1.32922 
1.32922 
1.32922 



2.727 
.527 



.33968 

.33853 
.33846 
.33845 
.33845 
.33845 



a 



2.806 
.005 



2.801 
.015 



2.803 
.020 



(oo, oo) 



1.5 



1.33333 .33333 



Table 5 

The mass /u(oo), dilaton charge D and Hawking Temperatur T/T5 for the S0(3) EYMD 
black hole solutions of the sequences + n) with < j < 2, 1 < n < 5 for 
dilaton coupling constant 7 = 1 and horizon xu — 1. For each sequence the limiting 
values /x4,(oo), -D^, and {T/Ts)^^ are given. The masses Unioo) are approximated by 
/Xoo(oo) — iin{oo) = a^e~"^" (and analogously D and T/Ts), with the constants a, a 
also shown. The last hne contains the values obtained for //^(oo), Di^ and {T/Ts)lo 
in the limit j — > 00. 
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// 


J = 


1 


2 




dfi 






1 


2.47869 


1.86167 


1.69430 


.094 


2.046 






2 


1.86074 


1.64685 


1.58248 


.040 


1.714 






3 


1.64914 


1.56969 


1.53269 


.024 


1.608 


7 = 


= 0.1 


4 


1.57468 


1.53460 


1.50527 


.019 


1.565 






5 


1.53986 


1.50966 


1.48215 


.018 


1.538 








.034 


.019 


.017 










dj 


1.821 


1.663 


1.611 










1 


2.29374 


1.34077 


.89474 


.404 


2.009 






2 


1.43827 


.94746 


.65374 


.371 


1.373 






3 


1.01694 


.69874 


.48400 


.367 


1.009 


7 = 


= 0.5 


4 


.75328 


.52643 


.36697 


.361 


.755 






5 


.56598 


.39571 


.27491 


.364 


.570 








.293 


.284 


.283 












2.443 


1.640 


1.133 










1 


2.18639 


.85307 


.34698 


.920 


2.171 






2 


1.02365 


.42121 


.17227 


.891 


1.025 






3 


.50393 


.20783 


.08377 


.897 


.506 


7 


= 1 


4 


.24955 


.10352 


.04215 


.889 


.250 






5 


.12364 


.05102 


.02055 


.897 


.124 








.716 


.704 


.706 












4.373 


1.722 


.704 










1 


2.62983 


.68046 


.16794 


1.376 


2.651 






2 


.91273 


.23353 


.05661 


1.390 


.921 






3 


.31115 


.07745 


.01818 


1.420 


.314 


7 


= 2 


4 


.10281 


.02540 


.00605 


1.416 


.103 






5 


.03357 


.00824 


.00193 


1.428 


.034 








1.091 


1.105 


1.117 












8.001 


2.096 


.520 







Table 6 

The location of the innermost node z^^ for the SO (3) EYMD regular solutions with 
node structure (j, j-\-n) with j <2 and n < 5 for the dilaton coupling constants 7 = 0.1, 
0.5, 1, 2. The nodes are approximated by z^\j) = dje~'^^^ and Zj^\n) = dne~'^"\ with 
the constants dj, Cj and dn, Cn also shown. 
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10 < xh < oo 





j = 
fiioo) D T/Ts 


J = l 
/i(cx)) D T/Ts 


J = 2 
/i(oo) D T/Ts 


6 

b 

a 


1.006 1.009 2.017 
2.553 5.302 5.439 
2.795 2.804 2.795 


1.008 1.016 2.024 
0.071 0.148 0.157 
2.807 2.804 2.808 


1.008 1.018 2.025 
0.002 0.003 0.004 
2.808 2.697 2.810 



Table 7 

The constants 5, 6, a of the least square fits for A/i^ 



/ioo(oo) -/i„(oo) 



(and analogous for D and T/Ts) for the S0(3) EYMD sequences (0, n) 
(2, 2 + n) for 7 = 1 and > 10. For < 1 see . 



(1, 1 + ra) and 
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atter Function (x 
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10 



Figure la: The SU(2) EYMD gauge field functions Wn{x) with node number n = 1 
(dotted), n = 3 (dashed), n = 5 (dot-dashed) and n = 7 (tripledot-dashed) for the 
black hole solutions with horizon x^i = 0.01 and dilaton coupling constant 7 = 1. The 
solid line shows the gauge field function of the limiting EMD solution with the same 
horizon and dilaton coupling constant. The thin vertical lines indicate the location of 
the innermost node of the corresponding regular solutions. 
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Dilaton Function ( 
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Figure lb: Same as Figure |T^for the dilaton functions (j)n{x). 
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trie Function (xj^=O.Ol) 
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Figure Ic: Same as Figure ^ for the metric functions Nn{x). 
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trie Function {x^=0 
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Figure Id: Same as Figure ^ for the metric functions An{x). 



46 



Charge (xh = 0.01) 
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Figure le: Same as Figure ^ for the squared charge functions -P^(a;). Also the function 
(1 + W7(x))^ is shown (thin tripledot-dashed). 
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alter Function ( 
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Figure 2a: The SU(2) EYMD gauge field functions Wn{x) with node number n — 1 
(dotted), n — 3 (dashed), n — 5 (dot-dashed) and n — 7 (tripledot-dashed) for the 
regular solutions with dilaton coupling constant 7 = 1. The solid line shows the gauge 
field function of the "extremal" EMD solution with the same dilaton coupling constant. 
The thin vertical lines indicate the location of the innermost nodes. 
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atter Function 
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Figure 3: The SU(2) EYMD gauge field functions W7{x) for tlie regular solution 
(dotted) and the black hole solutions with horizon xn = 0.01 (dashed), = 0.1 (dot- 
dashed) and xh — 1 (tripledot-dashed) for the dilaton coupling constant 7 = 1. The 
solid hne shows the gauge field function of the "extremal" EMD solution. 
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Figure 4a: The deviation at the horizon Awnixa) = \wn{x}i) — Woo{xu) \ for the SU(2) 
EYMD gauge field functions with up to five nodes as a function of the horizon xh for 
the dilaton coupUng constant 7 = 1. The vertical lines indicate the location of the 
innermost and second nodes of the corresponding regular solutions. The horizontal 
lines show the deviations of 0.75, 0.27, 0.05 and 0.008 from the limiting value. 
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Figure 4b: The relative deviation at the horizon A0„(a;H) = \{4>n{xu) — 
(f^ooix-n)) / (t>oo{x-ii)\ for the SU(2) EYMD dilaton functions with up to five nodes as 
a function of the horizon xh for the dilaton coupling constant 7=1. The vertical lines 
indicate the location of the innermost and second nodes of the corresponding regular 
solutions. 
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Hawking Temperature 




Figure 5: The inverse Hawking temperature /? as a function of the mass /i(oo) for the 
SU(2) EYMD solutions with node number n — 1 (dotted), n — 2 (tripledot-dashed) , 
n — 3 (dashed) and n = 4 (long dashed) for dilaton coupling constants 7 = 1 and 
7 = 0. The solid lines show the inverse Hawking temperature of the EMD and RN 
solutions, respectively. 
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Figure 6: The logarithm of the absolute deviation from the limiting solution A//„ = 
A*oo(oo) — A*n(oo) for the SU(2) EYMD masses as a function of the node number n with 
dilaton coupling constant 7 = 1 for the regular solution (solid line) and the black hole 
solutions (xh = 0.01, 0.1, 1, 10, 100 and 1000). The dotted lines indicate the least 
square fits with fit constants given in Table 3. 
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Matter Function (xh = 0.01) 
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Figure 7a: The SU(2) EYMD gauge field function tyy of the black hole solutions with 
horizon xu — 0.01 and dilaton couphng constants 7 = (dotted), 7 = 0.5 (dashed), 
7 = 1 (dot-dashed) and 7 = 2 (tripledot-dashed). The solid line shows the gauge 
field function of the EMD solution. The thin vertical lines indicate the location of the 
innermost node of the corresponding regular solutions. 
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Metric Function (xh=0.01) 




log(x) 



Figure 7b: The SU(2) EYMD metric functions Nrlx) of the black hole solutions with 
horizon = 0.01 and dilaton coupling constants 7 = (dotted), 7 = 0.5 (dashed), 7 = 
1 (dot-dashed) and 7 = 2 (tripledot-dashed) . The sohd hne shows the metric function 
of the Reissner-Nordstr0m solution and the long dashed line the metric function of the 
EMD solution with 7 = 0.5. For 7 = 1 and 2 the functions A^7(a;) and their limiting 
EMD functions fall on top of each other and are indistinguishable. The thin vertical 
lines indicate the location of the innermost node of the corresponding regular solutions. 
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Hawking Temperature 




Figure 8: The inverse Hawking temperature /? as a function of the mass f^{oo) for the 
SU(2) EYMD solutions with node number n — 1 (dotted), n — 2 (tripledot-dashed) , 
n — 3 (dashed) and n = 4 (long dashed) for the dilaton coupling constants 7 = 0.5 
and 7 = 2. The solid lines show the inverse Hawking temperature of the corresponding 
limiting EMD solutions. 
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Matter Function (5,5) 




log(x) 

Figure 9: The S0(3) EYMD gauge field functions K^{x) and H^{x) for the regular 
solution (dotted) and black hole solutions with horizon xu = 0.01 (dashed), xn = 0.02 
(dot-dashed) and — 0.03 (solid) with node structure (5, 5) and dilaton coupling 
constant 7 = 1. 
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Figure 10: The critical horizon a;g(7) as a function of the dilaton coupling constant 
7 of the S0(3) EYMD solutions with node structure (n, n) (dotted) together with the 
location of the innermost nodes of the corresponding regular solutions (dashed). Also 
shown are the locations of the innermost nodes of the scaled SU(2) solutions (dot- 
dashed) with the same node structure. 
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Matter Function (x 
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Figure 11a: The SO (3) EYMD gauge field functions Kn{x) and Hn{x) of the solutions 
with node structure (0,n) for dilaton coupling constant 7=1 and horizon = 0.2, 
for n = 1 (dotted), n = 3 (dashed) and n = 5 (dot-dashed). The solid line shows the 
gauge field function of the limiting EMD solution with magnetic charge — ?> and 
the same dilaton coupling constant and horizon. The thin vertical lines indicate the 
location of the innermost node of the corresponding regular solutions. 
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Metric Function {x^=0.2) 




EMD 



1 - 



n=l \ n=3 
\ 

\ 

\ 



(7 = 1) 



n = 5 



2 



2 

logfx 



4 



6 



Figure 11c: The same as Figure |lla| for the squared charge functions 
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Figure 12: The SO (3) EYMD gauge field functions K5{x) and H5{x) with node 
structure (0, 5) and dilaton couphng constant 7 = 1 for the regular solution (dotted) 
and black hole solutions = -5 (dashed), x-^ = 1 (dot-dashed) and xh = 2 (tripledot- 
dashed). The solid line shows the gauge field function of the "extremal" EMD solution 
with magnetic charge = 3 and the same dilaton coupling constant. 
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Hawking Temperature (0,n) 
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Figure 13: The inverse Hawking temperature /3 as a function of the mass f^{oo) for the 
S0(3) EYMD solutions with node structure (0, n) for the dilaton couphng constants 
7 = and 7 = 1 and node number n = 1 (dotted), n = 2 (tripledot-dashed) , n — 3 
(dashed) and n = 4 (long dashed). Also shown is the inverse Hawking temperature of 
the limiting EMD and RN solutions (solid). 
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Matter Function {x^ = 0.2) 




Figure 14a: The S0(3) EYMD gauge field functions Kn{x) and Hn{x) of the solutions 
with node structure (1,1 + n) for the dilaton coupling constant 7 = 1 and horizon 
xn = 0.2, for n = 1 (dotted), n = 3 (dashed) and n = 5 (dot-dashed). The solid line 
shows the gauge field function of the limiting charged SU(3) EYMD solution with one 
node, j — 1, and magnetic charge P'^ — 3. The thin vertical fines indicate the location 
of the innermost node of the corresponding regular solutions. 
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trie Function {x^=0.2) 
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Figure 14b: The same as Figure |14a| for the metric functions Nn{x). 
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Figure 14c: The same as Figure |14a| for the squared charge functions 
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Figure 15: The inverse Hawking temperature /? as a function of the mass ^{oo) for the 
S0(3) EYMD solutions with node structure (1, 1 + n) for the dilaton couphng constants 
7 = and 7 = 1 and node number n = 1 (dotted), n = 2 (tripledot-dashed) , n = 3 
(dashed) and n — 4 (long dashed). Also shown is the inverse Hawking temperature 
of the limiting charged SU(3) EYMD and EYM solutions with one node, j — I and 
magnetic charge P"^ — 3. 
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Matter Function {x^ = 0,2) 
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Figure 16a: The S0(3) EYMD gauge field functions Kn[x) and Hn{x) of the solutions 
with node structure {2,2 + n) for the dilaton coupling constant 7 = 1 and horizon 
xh = 0.2, for n = 1 (dotted), n = 3 (dashed) and n = 5 (dot-dashed). The solid line 
shows the gauge field function of the limiting charged SU(3) EYMD solution with two 
nodes, j = 2, and magnetic charge = 3. 
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Hawking Temperature 
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Figure 17: The inverse Hawking temperature /5 as a function of the mass fi{oo) for 
the charged SO (3) EYM and EYMD solutions with magnetic charge = 3 for dilaton 
couphng constants 7 = and 7 = 1 and node number j = I (dotted), j = 2 (tripledot- 
dashed), j = 3 (dashed) and j = 4 (long dashed). Also shown is the inverse Hawking 
temperature of the RN and EMD solutions with magnetic charge = 4 (sohd) and 
P2 = 3 (thin dotted). 
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